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1 Introduction 



Abstract 

in 

The latest results from the Soudan 2, MACRO and SuperKamiokande experiments on atmospheric 
■ neutrino oscillations are summarised and discussed. In particular a discussion is made on the Monte Carlo 
^ ' simulations used for the atmospheric neutrino flux. 

< 

<N 

. Atmospheric neutrinos are generated in the decays of hadrons produced in high energy cosmic ray (CR) 
interactions. A high energy primary cosmic ray, proton or nucleus, interacts in the upper atmosphere 
CN| ■ producing a large number of pions and kaons, which decay yielding muons and muon neutrinos; also the 
muons decay yielding muon and electron neutrinos. The ratio of the numbers of muon to electron neutrinos 
is about 2 and N v /N-& ~ 1. Atmospheric neutrinos are produced in a spherical surface at about 10-20 km 
above ground and they proceed towards the earth. 

If neutrinos have non-zero masses, one has to consider the weak flavour eigenstates v e , v^, v T and the 
mass eigenstates v\, i>2, v%. The weak flavour eigenstates Vi are linear combinations of the mass eigenstates 
v m through the elements of the mixing matrix Ui m : 

Oh! n = V, u lm v m (l) 

For 2 flavour eigenstates (V^, v T ) and 2 mass eigenstates (yi, ^3) 

^23 



^ . I = v-i cos 6>23 + v-i sin 9 2o 

\ v T = — v 2 sin 823 + vi cos 623 

where #23 is the mixing angle. The survival probability of a beam is 



(2) 



/1.27Am 2 • L\ 

P{y„ -» Vp) = 1 - sin 2 20 23 sin 2 f J (3) 

where Am 2 = m| — m\, L is the distance travelled by the v from production to detection. The probability 
for the initial to oscillate into a v T is 

(1 27Ato 2 • L\ 
^ J (4) 

Atmospheric neutrinos are well suited for the study of neutrino oscillations, since they have energies from 
a fraction of GeV up to more than 100 GeV and they travel distances L from few tens of km up to 13000 
km; thus LjE v ranges from ~ 1 km/GeV to 10 5 km/GeV. Moreover one may consider that there are two 
identical sources for a single detector: a near one (downgoing neutrinos) and a far one (upgoing neutrinos). 
Atmospheric neutrinos are particularly useful to study oscillations for small Am 2 , and matter effects can 
be studied with their high energy component. 
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Figure 1: Longitudinal cross section of the Soudan 2 detector and observed event topologies: (a) Fully 
Contained Events (E u ) ~ 1 GeV, (b) Partially Contained Events {E v ) ~ 6 GeV, (c) In-down muons 
(E u ) ~ 2.4 GeV, (d) Up-stopping muons (E„) ~ 6.2 GeV. 



The early water Cherenkov detectors and the tracking calorimeters measured and z/ e charged 
current interactions. The results were expressed in terms of the double ratio R' — Robs/ RmCi where 
Robs = {Ny^/N^oks is the ratio of observed fj, and e events and Rmc — (^/^Imc is the same ratio 
for Monte Carlo (MC) events. The R 1 double ratios from 1MB [1] and Kamiokande [2] were smaller than 
expectations, while the NUSEX [3] and Frejus [4] R' agreed with expectations. The Baksan [5] scintillation 
telescope detected upthroughgoing muons at the expected rate but gave indications of an anomalous angular 
distribution. Later, the Soudan 2 tracking and shower calorimeter detector confirmed the anomaly in the 
Vft/vg double ratio for contained events [6]. MACRO reported in 1995 a measurement of upthroughgoing 
muons coming from of (E u ) ~ 50 GeV, in which there was an anomalous zenith distribution and a deficit 
in the total number of observed upgoing muons [7]. SuperKamiokande (SK) confirmed the anomalous 
double ratio and provided a wealth of informations for sub-GeV and multi-GeV v^, v e and for higher energy 
upthroughgoing muons and stopping muons. In 1998 Soudan 2, MACRO and SK provided strong indications 
in favour of < — > v T oscillations [8-11]. After 1998 new results were presented by the 3 experiments [12-14]. 
Here we shall review their results. 



2 Atmospheric neutrino flux calculations 

In the past use was made of unidimensional Monte Carlo codes, Bartol96 [15] and HKKM95 [16]. Recently 
new improved MC predictions for neutrino fluxes were made available by the HKKM01 [17] and FLUKA 
[18] groups. They include three dimensional calculations of hadron production and decays and of neutrino 
interactions, improved hadronic model and new fits of the primary CR flux. The two MCs yield predictions 
for the non oscillated and oscillated fluxes equal to within few %. The shapes of the angular distributions 
for oscillated and non oscillated Bartol96, HKKM95, new FLUKA and new HKKM01 fluxes are the same 
to within few %. The absolute values of the MACRO upthroughgoing muon data are about 25% higher 
than those predicted by the new FLUKA and HKKM01 MC [13], while the shapes of the oscillated and non 
oscillated angular distributions differ by no more than 5%, see Sec. 4. 

A similar situation is found in the SK data [14], see Sec. 5. The clcctron-likc events were in agreement 
with the HKKM95 [16] MC predictions in absence of oscillations, while they are higher than the HKKM01 
[17] non oscillated MC, Fig. 7. For the muon-like events, the new MC predictions are low for the SK data, 
especially for the high energy upthroughgoing events [14]. Previous comparisons between the SK muon data 
and the HKKM95 [16] predictions showed a global deficit of events and a zenith distribution in agreement 
with < — ► v T oscillations [14]. 

The difference between the new and old MC predictions is very probably due to the use of a new fit of 
the cosmic ray data [19]. 

Recent results by the L3C and BESS experiments [20] on the primary cosmic ray fit show good agreement 
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Figure 2: Zenith angle distribution for Hi-Res e— flavour (a) and fi— flavour (b) events. The points with 
error bars are the data, the dashed and solid histograms are the predicted non oscillated and oscillated 
v distributions, respectively, the dotted histograms are the contribution of the rock background, (c) 
Distribution in log LjE v for the Hi- Res Fully Contained events (black crosses) compared with the MC 
predictions for oscillations (solid histogram), no oscillations (dashed histogram) and saturated oscillations 
(dotted histogram). 



with the Bartol96 and HKKM95 predictions and a disagreement with the new fit of the cosmic ray data 
[19]. 

The calculations of HKKM01 and FLUKA are in good agreement when using the old fit or the new fit 
to the primary CR flux; this confirms the improvement in the hadronic model. 

In SK the evidence for neutrino oscillations lies in the shapes of the angular distribution and in the ratio 
of /z/e data. In MACRO the evidence is due mainly to the shape of the high energy angular distribution 
and the prediction is the same in all simulations. 

The MACRO data suggest that the FLUKA normalisation should be raised by ~ 25% at E v ~ 50 GeV 
and by 12% at E v ~ 2 — 3 GeV [13]. Similar conclusions are reached by SK [14]. 

3 Results from the Soudan 2 experiment 

The Soudan 2 experiment used a modular fine grained tracking and showering calorimeter of 963 t, located 
2100 m.w.e. underground in the Soudan Gold mine in Minnesota. The bulk of the mass consisted of 
1.6 mm thick corrugated steel sheets interleaved with drift tubes. The detector was surrounded by an 
anticoincidence shield. Figure 1 shows a longitudinal cross section of the apparatus and the topologies of 
the events observed. The final analysis used the Fully Contained High- Resolution events [12], Table 3. 





Data 


MCno osc [15] 


Track 


101.9 ± 12.7 


193.1 


Showers 


146.7 ± 12.5 


179.0 



Table 1: Soudan 2 Hi- Res data. 



The fully contained events consist mostly of quasi-elastic neutrino reactions, but include a background 
of photons and neutrons from cosmic ray muon interactions in the surrounding rock. The track and shower 
events for a 5.9 kt-yr exposure are summarised in Table 3, where they are compared with MC predictions 
based on the Bartol96 neutrino flux [15]. 

After corrections for background and selecting a high resolution (Hi-Res) sample of events, the Soudan 
2 double ratio for the whole zenith angle range (—1 < cosO < 1) is R' = {N IJL /N^)DATA/{N l i/N e )MC = 
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Figure 3: Cross section of the MACRO detector. Event topologies induced by interactions in or around 
the detector. IU M = semicontained Internal Upgoing fj,; ID ;t = Internal Downgoing /i; UGS M = Upgoing 
Stopping fi; Upthroughgoing = upward throughgoing fi. 



0.69 ± 0.12, consistent with muon neutrino oscillations. 

The v e data agree with the no oscillation MC predictions, while the data are lower, except in the 
forward zenith bin. The double peak structure arises from the acceptance of the apparatus. The roughly 
interpolated 90% C.L. allowed region in the sin 2 29 — Am 2 plane, computed using the Feldman-Cousins 
method [21] is shown in Fig. 8, where it is compared with the allowed regions obtained by the SK and 
MACRO experiments. 

4 Results from the MACRO experiment 

MACRO was a large area multipurpose underground detector designed to search for rare events and rare 
phenomena in the penetrating cosmic radiation. It was located in Hall B of the Gran Sasso Lab at an 
average rock overburden of 3700 m.w.e.; it started data taking with part of the apparatus in 1989; it was 
completed in 1995 and was running in its final configuration until the end of 2000. The detector had global 
dimensions of 76.6 x 12 x 9.3 m 3 ; vertically it was divided into a lower part, which contained 10 horizontal 
layers of streamer tubes, 7 of rock absorbers and 2 layers of liquid scintillators, and an upper part which 
contained the electronics and was covered by 1 scintillator layer and 4 layers of streamer tubes. The sides 
were covered with 1 vertical scintillator layer and 6 of limited streamer tubes [22]. 

MACRO detected upgoing v^s via charged current interactions — > /z; upgoing muons were identified 
with the streamer tube system (for tracking) and the scintillator system (for time-of- flight measurement). 
The events measured and expected for the 3 measured topologies, Table 4, and the L/E v distribution, 
Fig. 5a, deviate from MC expectations without oscillations; the deviations point to the same < — ► v T 
oscillation scenario [13]. 

Upthroughgoing muons (E^ > 1 GeV) come from interactions in the rock below the detector of 
with (E v ) ~ 50 GeV. The MC uncertainties arising from the neutrino flux, cross section and muon 

propagation on the expected flux of muons are estimated to be ~ 17%; this systematic error on the 

upthroughgoing muons flux is mainly a scale uncertainty. 

In order to verify that different flux simulations affect the zenith distribution at the level of only a few 

percent (while there is an effect of the order of ~ 25% on the event rates) MACRO compared the predictions 

of the Bartol96 [15], FLUKA [18] and HKKM01 [17] MCs. In Fig. 4a the MACRO data are compared with 
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Events 


MCno osc [15] 


R = Data/MCno osc 


Upthr. 


857 


1169 


0.73 


IU 


157 


285 


0.55 


ID+UGS 


262 


375 


0.70 



Table 2: MACRO events. 




Figure 4: (a) Comparison between the zenith distribution of the MACRO upthroughgoing muons and the 
oscillated MC predictions given by Bartol96 (solid curve), HKKM01 (dash-dotted line), FLUKA fitted to the 
new CR measurements (dashed curve) and FLUKA with the old CR fit (dotted curve) . Zenith distributions 
for (b) IU and for (c) ID+UGS MACRO events (black points) compared with the no oscillation Bartol96 
MC (dashed line with a scale error band) and with the < — ► v T predictions with Am 2 = 2.3 • 10~ 3 eV 2 
and maximal mixing. 



the oscillated Bartol96, the new HKKM01 and FLUKA calculations using the new CR fit. All predicted 
curves are for maximal mixing and Am 2 = 0.0023 eV 2 . 

Low energy events. Semicontained upgoing muons (IU) come from interactions inside the lower 
apparatus. The average parent neutrino energy for these events is ~ 2 — 3 GeV. Up stopping muons (UGS) 
are due to external interactions yielding upgoing muons stopping in the detector; the semicontained 
downgoing muons (ID) are due to downgoing z/^'s with interaction vertices in the lower detector. The lack 
of time information prevents to distinguish between the two subsamples. An almost equal number of UGS 
and ID events is expected. The average parent neutrino energy for these events is ~ 2 — 3 GeV. The number 
of events and the angular distributions are compared with the MC predictions without oscillations in Table 
4 and Fig. 4b, c. The low energy data show a uniform deficit of the measured number of events over the 
whole angular distribution with respect to the Bartol96 predictions. 

< ► v T against < ► v sterile- Matter effects due to the difference between the weak interaction 

effective potential for muon neutrinos with respect to sterile neutrinos, which have null potential, yield 
different total number and different zenith distributions of upthroughgoing muons. The measured ratio 
between the events with —1 < cosO < —0.7 and with —0.4 < cosO < was used [13]. In this ratio most 
of the theoretical uncertainties on neutrino flux and cross sections cancel. Combining the experimental and 
theoretical errors in quadrature, a global uncertainty of 6% is obtained. The measured ratio is R m eas = 1-38, 
to be compared with R T = 1.61 and R s teriie = 2.03. One concludes that < — > v s teriie oscillations (with 
any mixing) are excluded at the 99.8% C.L. compared to the < — ► v T channel with maximal mixing and 
Am 2 = 2.3 • 10" 3 eV 2 . 
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Figure 5: (a) Ratio Data/MCno osc as a function of the estimated L/E v for the upthroughgoing muon 
sample (black points). The solid line is the MC expectation assuming Am 2 = 2.3 • 1CP 3 eV 2 and sin 2 29 = 1. 
The last point (empty circle) is obtained from the IU sample, (b) Interpolated 90% C.L. contour plots 
of the allowed regions in the sin 2 29 — Am 2 plane for the MACRO data using only the ratios i?i , i?2 , R3 
(continuous line) and adding also the information on the absolute values i?4, -R5 (dotted line). 



energy estimate by Multiple Coulomb Scattering of muons. Since MACRO was not equipped 
with a magnet, the only way to estimate the muon energy is through their Multiple Coulomb Scattering 
(MCS) in the absorbers. Two analyses were performed [23]. The first was made studying the deflection of 
upthroughgoing muons using the streamer tubes in digital mode. This method had a spatial resolution of ~ 1 
cm. The second analysis was performed using the streamer tubes in "drift mode" [23] . The space resolution 
was ~ 3 mm. For each muon, 7 MCS variables were defined and given in input to a Neural Network, 
previously trained with MC events of known energy crossing the detector at different zenith angles. The 
output of this program gave the muon energy estimate event by event. The sample of upthroughgoing 
muons was separated in 4 subsamples with average energies of 12, 20, 50 and 100 GeV. The ratios 
Data/MCno osc as a function of \og 1Q {L/E v ) obtained from upthroughgoing muons are plotted in Fig. 5a; 
they are in agreement with the < — ► v T oscillation hypothesis [13]. 

New determination of the oscillation parameters. In previous analyses MACRO fitted the shape 
of the upthroughgoing muon zenith distribution and the absolute flux compared to Bartol96. This yielded 
Am 2 = 2.5 ■ 10 -3 eV 2 and maximal mixing [13]. Later, in order to reduce the effects of systematic 
uncertainties in the MC simulations, MACRO used the following three independent ratios. It was checked 
that FLUKA, HKKM01 and Bartol96 Monte Carlo simulations yield the same predictions to within ~ 5%. 

(i) High Energy Data: zenith distribution ratio: R\ — N vert /Nh or (ii) High Energy Data, v energy 
measurement ratio: R2 — Ni ow /Nhi g h (hi) Low Energy Data: R3 = (Data/MC)iu /(Data/MC)m+UGS- 
The no oscillation hypothesis has a probability P ~ 3 ■ 10~ 7 and is thus ruled out by ~ 5o\ By fitting the 3 
ratios to the < — > v T oscillation formulae, MACRO obtained sin 2 29 = 1, Am 2 = 2.3 • 10~ 3 eV 2 and the 
allowed region indicated by the solid line in Fig. 5b. There is a good consistency between the old and new 
methods. 

Using the Bartol96 flux, it is possible to add the information on the absolute flux values of the 
(iv) High energy data (systematic error ~ 17%) i?4 = N meas /N MC - (v) Low energy semicontained 
muons (scale error 21%) R 5 = N meas /N MC - 

These informations reduce the area of the allowed region in the sin 2 29 — Am 2 plane, as indicated by the 
dashed line in Fig. 5b. The final MACRO best fit is Am 2 = 2.3 • 10~ 3 eV 2 and sin 2 26 = 1 (6cr significance). 
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Figure 6: Sharpness of Chcrcnkov rings produced in the SK detector by (a) an electron and (b) a muon. 





Data 


MCno osc [17] 


Sub-GeV 1-ring 


3227 


4213 


Sub-GeV multi-ring 


208 


323 


Multi-GcV 1-ring 


651 


900 


Multi-GeV multi-ring 


439 


712 


PC n 


647 


1034 


Up-stop 


418 


721 


Upthr. 


1842 


1684 



Table 3: SK fj, data taken from 1996 to 2001. 



5 Results from the SuperKamiokande experiment 

SuperKamiokande [14] is a large cylindrical water Cherenkov detector of 39 m diameter and 41 m height 
containing 50 kt of water (the fiducial mass of the detector for atmospheric neutrino analyses is 22.5 kt); it 
was seen by 11146, 50-cm-diameter inner-facing phototubes. The 2 m thick outer layer of water, acting as 
an anticoincidence, was seen by 1885 smaller outward-facing photomultiplicrs. The ultra pure water has a 
light attenuation of almost 100 m. The detector is located in the Kamioka mine, Japan, under 2700 m.w.e. 
SK took data in its full configuration from April 1996 till November 2001, when an accident happened. It 
resumed operation with about half of PMTs in October 2002. 

Atmospheric neutrinos are detected in SK by measuring the Cherenkov light generated by the charged 
particles produced in the neutrino CC interactions with the protons and oxygen nuclei. The large detector 
mass and the possibility of clearly defining a large inner volume allow to collect a high statistics sample 
of fully contained events (FC) up to relatively high energies (up to ~ 5 GeV). The FC events have both 
the neutrino vertex and the resulting particle tracks entirely within the fiducial volume; they yield rings 
of Cherenkov light on the PMTs. Fully contained events can be further subdivided into two subsets, the 
so-called sub-GeV and multi-GeV events, with energies below and above 1.33 GeV, respectively. 
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Another sub-sample, defined as the partially contained events (PC), is represented by those CC 
interactions where the vertex is still within the hducial volume, but at least a primary charged particle, 
typically the muon, exits the detector without releasing all of its energy. In this case the light pattern is a 
filled circle. For these events the energy resolution is worse than for FC interactions. Upward- going muons 
(UPMU), produced by neutrinos coming from below and interacting in the rock, are further subdivided into 
stopping muons {{E v ) ~ 7 GeV) and upthroughgoing muons ({E u ) ~ 70 -f- 80 GeV), according to whether 
or not they stop in the detector. 
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Figure 7: SK data taken from 1996 till 2001 with the detector in full configuration (1489 days for FC+PC 
events and 1646 days for upgoing fi). Zenith distributions for e-like and /i-like sub-GeV and multi-GeV 
events, for partially contained events and for upthroughgoing and stopping muons (black points). The 
boxes are the no oscillation HKKM01 predictions, the solid lines refer to < — > v r oscillations with 
maximal mixing and Am 2 = 2.4 • 10 -3 eV 2 . 



Particle identification in SuperKamiokande is performed using likelihood functions to parametrise the 
sharpness of the Cherenkov rings, which are more diffuse for electrons than for muons, Fig. 6. The 
algorithms are able to discriminate the two flavours with high purity (of the order of 98% for single track 
events). The zenith angle distributions for e-like and /x-like sub-GeV and multi-GeV events, for PC events 
and for upward throughgoing or stopping muons are shown in Fig. 7. These data were taken from 1996 
till 2001 with the detector in full configuration (1489 days for FC+PC events and 1646 days for upgoing 
/i). The data and MC behaviour shows the problem with the new HKKM01 MC discussed in Sec. 2. The 
number of measured and expected /i-like events are summarised in Table 5. The new data, taken in 2003 
with about half of the PMTs and referring to 311 days for FC+PC events and 243 days for upgoing muons, 
show the same behaviour [14] as the older data. 

New analyses have been performed leaving free the normalisation. The last value for the double ratio 
B! reported by SK is 0.658 ± 0.016 stat ± 0.032 sys for the sub-GeV sample and 0.702 ± 0.031 stQt ± 0.099 sys 
for the multi-GeV sample (both FC and PC). 

SK used also a selected sample of events with good resolution in LJE V , to search for the dip in the 
oscillation probability expected when the argument of the second sine-squared term in Eq. (4) is tt/2. A dip 
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Figure 8: (a) Ratio of the data to the MC events without neutrino oscillation (black points) as a function 
of the reconstructed L/E v compared with the best-fit expectation for 2-flavour ia, < — > v T oscillations (solid 
line). The error bars are statistical only, (b) 90% C.L. allowed regions for Vu < — ► v T oscillations obtained 
by the SK, MACRO and Soudan 2 experiments. 



in the L/E u distribution is observed at L/E v ~ 500 km/GeV, see Fig. 8a. This is another proof in favour 
of v oscillations and a further constraint on Am 2 . Alternative models that could explain the zenith angle 
and energy dependent deficit of the atmospheric muon neutrinos are disfavoured, since they do not predict 
any dip in the L / E v distribution [14] . 

Interpreting the /x-like event deficit as the result of < — ► v r oscillations in the two-flavour mixing 
scheme, SK computed an allowed domain for the oscillation parameters [14], see Fig. 8b. The events 
were binned in a multi-dimensional space defined by particle type, energy and zenith angle, plus a set of 
parameters to account for systematic uncertainties. The best fit using FC, PC, UPMU and MRING events 
[14] corresponds to maximal mixing and Am 2 = 2.4 • 10~ 3 eV 2 . In Fig. 8b, the 90% SK allowed region 
in the sin 2 29 — Am 2 plane is compared with the MACRO and Soudan 2 ones. The limit lines represent 
smoothed interpolations and are qualitative. 



6 Conclusions 

The atmospheric neutrino data strongly favour « — > v T oscillations with maximal mixing and Am 2 = 
0.0023 - 0.0052 eV 2 (Soudan 2: 0.0052, MACRO: 0.0023, SK: 0.0024 eV 2 ). MACRO and SK exclude 
Up < — ► v st eriie oscillations at the level of 99.9%; SK excludes also < — ► v e oscillations and some exotic 
processes [14, 24]. The K2K [25] long baseline experiment, using the beam produced at KEK and 
detected by SK, supports the < — > v T oscillations with Am 2 = 0.0027 eV 2 . SK finds the predicted dip in 
the L/E u distribution. 

It has been hypothesized that, besides the dominant mass neutrino oscillations, there could be sub- 
dominant oscillations due to possible Lorentz Invariance Violation (LIV) [26] (or violation of the equivalence 
principle). In this case, one would consider a mixing of flavour eigenstates and velocity eigenstates and 
estimate upper limits on the LIV parameters Sv/2 — (1)3 — v-i)/2 and sin 2 29 v . Preliminary analyses of SK 
and MACRO high energy muon data yield upper limits for Sv/2 at the level of ~ 10 -23 for low mixing 
angles and ~ 10~ 26 for large mixing angles [27]. 

Several other possibilities have been discussed, also possible v decays [27, 28]. 
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